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Abstract. A first-order phase transition is found in two types of intrinsic curvature models defined on 
dynamically triangulated surfaces of disk topology. The intrinsic curvature energy is included in the Hamil- 
tonian. The smooth phase is separated from a non-smooth phase by the transition. The crumpled phase, 
which is different from the non-smooth phase, also appears at sufficiently small curvature coefficient a. 
The phase structure of the model on the disk is identical to that of the spherical surface model, which was 
investigated by us and reported previously. Thus, we found that the phase structure of the fluid surface 
model with intrinsic curvature is independent of whether the surface is closed or open. 



PACS. 64.60.-i General studies of phase transitions - 68.60.-p Physical properties of thin films, nonelec- 
tronic - 87.16.Dg Membranes, bilayers, and vesicles 



1 Introduction 



Triangulated surfaces are one of the basic models to inves- 
tigate the physics for biological membranes and strings 1 
2,3,4|5]. Theoretical investigations, on the basis of renor- 
malization group strategy, have been made to reveal how 
the fluctuations of surfaces stiffen/soften the surface [B][7| 
8 , 9 J"0] . A considerable number of numerical studies have 
also been performed for understanding the phase struc- 
ture by using the surface simulations on fluid (or dy- 

"ll"] and on tethered (or fixed connectivity) surfaces [T§ll2Til 

model is defined on dynamically triangulated surfaces in 
the framework of the surface simulations, whereas the 
tethered model is defined on fixed connectivity surfaces. 
The surface simulations have relatively long history rather 
than that of another numerical technique Molecular Dy- 
namics simulations for the lipid models for biological mem- 
branes 37 , 38 39 , 40 . The extrinsic curvature energy is al- 
ways included in the surface Hamiltonian, because it is 
considered to smoother the surface. However, we must re- 
call that the intrinsic curvature energy such as the deficit 
angle term can also be assumed [4111421143*] in the Hamil- 
tonian because it can smoother the surface. 

It has been shown recently by Monte Carlo (MC) sim- 
ulations that the extrinsic curvature model of Helfrich 
and Polyakov-Kleinert 44,45,46 undergoes a discontin- 
uous transition between the smooth phase and the crum- 
pled phase on fixed connectivity spherical surfaces "2*8TI!M 
I3(J| . Moreover, the phase transition occurs independent of 



whether the Gaussian bond potential is included in the 
Hamiltonian or not "2""]. It was also shown that the tran- 
sition occurs independent of how the extrinsic curvature 
term is discretized [3T]. It has also been reported that 
a tethered surface model undergoes a first-order transi- 
tion, where the Hamiltonian includes an intrinsic curva- 
ture term [191121)112 ]"| . The transition in this model occurs 
independent of whether the surface is closed or open. In 
fact, the transition can be observed on a sphere [T""f|, on a 
disk [201 , and on a torus [21] • 

However, it seems that many points remain to be stud- 
ied in the surface models. We consider that one interesting 
dproblem is how the fluidity of lateral diffusion influences 
the phase transition in the model with some curvature 

agram of the intrinsic curvature model was studied on 
dynamically triangulated spheres [19|. however, it has not 
been studied yet on dynamically triangulated surfaces of 
disk topology. 

If a membrane model undergoes a phase transition on 
dynamically triangulated surfaces, we must have reasons 
why such transition occurs. In fact, it is widely accepted 
that we can see no phase transition in the fluid mem- 
branes. The standard argument for fluid membranes is 
that there is no long-range order on the surface at finite 
bending rigidity |4Tj . 

First of all, we must not forget that the transition 
seems to occur on relatively small surfaces |11II12II13"1|14| . 
Even if the model in this paper undergoes a transition 
on small surfaces, the transition may disappear on suf- 
ficiently large surfaces. Secondly, it is possible that the 
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transition occurs without the smooth phase: the transi- 
tion is not the one between the smooth phase and the 
crumpled phase but the one between the crumpled phase 
and the folded phase, which are respectively characterized 
by the Hausdorff dimension 2 < H < 3 and H > 3 in the 
case of self-intersecting surfaces. We have some other pos- 
sibilities of the origin including the argument of Helfrich 
that the thermal fluctuation increases the effective bend- 
ing rigidity 0H| . It is possible that the curvature coefficient 
is increased by the fluctuation of surfaces just like the case 
for the coefficient of the Gaussian curvature |35| , which is 
an intrinsic curvature as well as the deficit angle term in 
this paper. 

In this paper two types of fluid surface models are stud- 
ied by MC simulations on a disk with intrinsic curvature. 
The Hamiltonian of the first model is a linear combina- 
tion of the Gaussian bond potential and an intrinsic curva- 
ture term. The second model is a tensionless model, whose 
Hamiltonian is given by a linear combination of the intrin- 
sic curvature term and a hard-wall potential, which gives 
only an upper bound on the bond length. 

The simulations are performed under fixed total num- 
ber of vertices on the boundary of disk as well as on the 
interior of disk. As a consequence, the mean boundary 
length is fixed to be ub(1), where ub is the total num- 
ber of bonds on the boundary, and (£) is the mean bond 
length. 

It will be shown that both models undergo a first- 
order transition between the smooth phase and a non- 
smooth phase on relatively large surfaces. Thus, we can 
exclude the above-mentioned first and the second possibil- 
ities of the reason for the transition in the fluid membrane 
model with intrinsic curvature. As a consequence, the re- 
sult in this paper suggests that the fluctuation of surfaces 
does not always soften the surface. Thus we understand 
that fluidity of dynamically triangulated surfaces does not 
eliminate the long-range order whenever the surface model 
is defined with intrinsic curvature. This is in sharp con- 
trast to the case in the model with extrinsic curvature, 
where the surface is softened by the surface fluctuations. 

Combining the results in this paper and those in Ref. 
|19|. we will confirm that the transition is independent 
of whether the surfaces is closed or open in fluid and 
intrinsic-curvature models. We consider that this result is 
remarkable because phase transitions are always severely 
influenced by boundary conditions that restrict the value 
of variables in statistical models. 



2 Models 

The starting configuration for MC simulations is a flat and 
regularly triangulated disk, which is obtained by dividing 
the hexagon. Figure is a lattice of size AT = 61, which is 
obtained by dividing every edge of the original hexagon 
into 4-pieces. By dividing the edges into L-pieces we have 
a lattice of (N,N E ,N P ) = (3L 2 + 3L + 1,9L 2 + 3L,6L 2 ), 
where N,N E ,Np are the total number of vertices, the 
total number of edges, and the total number of plaquettes 
(or triangles). 




Fig. 1. A triangulated disk of size (N, N B , N P ) = (61, 156, 96), 
where N, Ne, Np are the total number of vertices, the total 
number of edges, the total number of plaquettes (or triangles) . 

The Gaussian bond potential Si, the intrinsic curva- 
ture S3, and the hard- wall potential V are given by 

5i = E ( x > - x i) 2 > 53 = E (*< - s o) 2 > 

v = Y,v{\x i -x j \), (1) 

where Y^Uj) ^ s the sum over all bonds (ij) connecting the 
vertices Xi and Xj, and Si in S3 is the sum of angles of 
the triangles meeting at the vertex i. We note that Yluj) 
in Si and V include the sum over the boundary bonds. 
As a consequence, the length of the boundary bonds can 
vary in the simulations. The value of the symbol <5o in S3 
is given by 

^ _ ( 2tt (internal vertices); 

1 7r(or 27t/3) (boundary vertices), ^ ' 

where 27r/3 is assigned to the six vertices, and tt to the 
remaining vertices on the boundary. Because of the def- 
inition of So in J21, we can call S3 in Eq.JTJ as a deficit 
angle term. 

The symbol V(\Xi—Xj\) in the hard-wall potential V 
in Eq. 0) is the potential between the vertices i and j, 
and is defined by 

y{\X t Xj\) — / ^ (otherwise). ^ 

The value of tq in the right hand side of Eq. © is fixed to 
r 2 = 1.08. Then we have (E( x i-Xj) 2 )/N ~ 3/2, which is 
satisfied when the Gaussian bond potential Si = ^ (Xi — 
Xj) 2 is included in the Hamiltonian, which includes no 
hard- wall potential such as V in Eq.Q. 

The constraint Xj \ <ro is necessary for the bond 
length to have a well-defined value if the Gaussian bond 
potential Si is not included in the Hamiltonian. In fact, 
the size of the surface grows larger and larger in the MC 
simulations if it were not for the constraint \Xi — Xj \ <rg 
in the hard-wall potential. As a consequence, the model 
seems to depend on a hidden length scale introduced by 
ro, and then the model appears to be ill-defined. However, 



H.Koibuchi: Phase structure of fluid membrane models with intrinsic curvature 



3 



we checked in that there is no ro-dependence on the 
results. This is a consequence of the scale invariant prop- 
erty of the model. Therefore, we use r% = 1.08 in the MC 
simulations. 

The partition function is defined by 

N 

Z{a) = J2 / II dX * ex P T )l ' ( 4 ) 

where N is the total number of vertices as described above, 
and denotes the summation over all possible triangu- 
lations T under a constraint which will be described in 
the next section. The expression S(X, T) shows that S 
explicitly depends on the variable X and T. The coeffi- 
cient a is the elastic modulus. The surfaces are allowed to 
self-intersect. The center of surface is fixed to remove the 
translational zero-mode. 

Two types of models, which we call model 1 and model 
2, are defined by Hamiltonians 

S(X, T) = Si+ aS 3 (model 1), (5) 
S(X, T) = V + aS 3 (model 2). (6) 

It should be noted that the Gaussian bond potential 
Si per vertex at the internal vertices differs from the one 
at the boundary vertices. In fact, the co-ordination num- 
ber is q — 6 at the internal vertices and q = 3 or q = 4 
at the boundary vertices. The intrinsic curvature S3 at 
the internal vertices can also be different from those at 
the boundary ones due to the same reason for Si. This 
difference of Si per vertex (or S3 per vertex) is typical 
to the model on the open disk, and such large difference 
cannot be seen in the model on closed surfaces such as a 
sphere or a torus. It must also be noted that the bound- 
ary is not completely free but shares both of the Gaussian 
bond potential (or the hard-wall potential) and the intrin- 
sic curvature energy in the models of this paper. 

We comment on a relation of S3 in Eq. Q to an intrin- 
sic curvature energy S3 = — J2i l°g(^/27r) in [T5|. This 
S3 in ^5] is closely related with the co-ordination de- 
pendent term S3 = — J2i l°g(<7i/6), which comes from the 
integration measure J\ i dXiqf 05] in the partition func- 
tion for the model on a sphere. S3 = — J2i l°g(<5i/27r) has 
the minimum value when Si = 2tt is satisfied for all i, 
and hence it becomes smaller on a smooth sphere than 
on a crumpled one. However, the flat configuration of 
the hexagonal lattice shown in Fig^ does not minimize 
S3 = — J2i l°g( ( 5i/27r), because it includes vertices such 
that Si 2ir. This is the reason why we use S3 in Eq.lJJJ, 
which becomes zero on such flat configuration as shown in 
Figd 

3 Monte Carlo technique 

The canonical Monte Carlo technique is used to update 
the variables X to a new X' so that X' = X + SX , where 
the small change SX is made at random in a small sphere 
111 R 3 . The new position X' is accepted with the probabil- 
ity Min[l,exp(-Z\S)], where AS = S(new) - S(old). The 



radius Sr of the small sphere is chosen at the beginning of 
the simulations to maintain the rate of acceptance rx for 
the X-update as 0.4 < r x < 0.6. 

The summation over T in Z of Eq. (jJJ is performed by 
bond flips. The bonds are labeled by sequential numbers. 
A bond is randomly chosen, and the flip is accepted with 
the probability Min[l, exp(— AS)]. The bond flip is applied 
only to internal bonds except the bonds whose both end 
points touch the boundary after the flip. The boundary 
bonds remain unflipped. 

iV-updates for X and TV-updates for T are consecu- 
tively performed, and these make one MCS (Monte Carlo 
Sweep). We use a random number called Mersenne Twister 
[5U| in the MC simulations. The above mentioned Metropo- 
lis procedure is applied only for X and T, which satisfy 
the constraint \Xi— Xj \ <ro in Eq.© for model 2. 

We use surfaces of size N = 817, N = 1657, N = 3367, 
and iV = 5167, which correspond to the partitions L = 16, 
L = 23, L = 33, and L = 41 respectively. The size of the 
surfaces is relatively larger than that of spherical surfaces 
used in [T5|. 

The total number of MCS after the thermalization is 
about 1.0 x 10 8 ~ 3.0 x 10 8 in the smooth phase close to 
the transition point of the surfaces of size N = 5167, N = 
3367, and TV =1657, and 1.6 xlO 8 for iV = 817. Relatively 
smaller number of MCS (1.0x10 s ~ 1.6xl0 8 ) (including the 
thermalization MCS) is done at non-transition point. The 
thermalization MCS is 0.5xl0 7 ~ 1.0x10 s , which depends 
on N and a. Very large number of thermalization MCS 
is sometimes needed close to the transition point, because 
there are smooth surfaces that start to collapse into non- 
smooth surfaces after a long MC run. Physical quantities 
are calculated at every 500 MCS after the thermalization. 



4 Results 

Snapshots of model 1 surfaces are shown in FigsEfa),EJl)), 
and|2fc), which are obtained at a = 2, a = 22, and a = 
28, respectively. Figure Eld) is a set of sections of the 
surface in FigEJc). The figures are shown in the same 
scale, and the size is N — 5167. It is easy to see that the 
surface in FigEJa) is crumpled (or folded), and the surface 
in FigEfb) is a branched polymer like one. The surface in 
FigGI c ) is obviously smooth, and it looks like a sphere 
with a hole. The surface in FigEJc) is separated from that 
in Fig[2b) by a first-order transition as we will see below. 

Snapshots of model 2 surfaces are shown in Figs Efa) > 
Gib), Etc), andEtd), which are obtained at a = 10, a = 32, 
a = 36, and a = 44, respectively. The figures are shown in 
the same scale, and the size is N = 5167. It is also easy 
to see that the surface in FigEIa) is almost folded, and 
that the surface in FigEfb) is a branched polymer like 
one, which is rather oblong than the one in FigEtb) of 
model 1. The snapshot in FigEIc) shows an intermediate 
and transient state between the non-smooth phase and the 
smooth phase. The reason why the surface in FigEIc) is 
considered to be an intermediate one is because the bend- 
ing energy S2 is relatively large, which will be shown in 
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(a) Model 1 (a = 2) 


(b) Model 1 (a = 22) 
















-•■X - 





(c) Model 1 (a = 28) (d) Sections of the surface in (c) 



Fig. 2. Snapshots of model 1 surfaces of size iV = 5167 at (a) 
a — 2 (crumpled), (b) a = 22 (non-smooth), and (c) a — 28 
(smooth), and (d) sections of the surface in (c). Figures are 
shown in the same scale. 

the next paragraph. We can see from Fig^c) that non- 
smooth part is almost shrunk. It is expected that the sur- 
face becomes a smooth surface after a long MC run. Be- 
cause the smooth phase is separated from the non-smooth 
phase by first-order transition, we expect that the configu- 
ration in Figl^c) eventually belongs to the smooth phase. 
On the contrary, the surface in Fig|3Jd) is almost smooth, 
and the non-smooth part seems to be shrunk. We note 
that the surface seems to be closed in the smooth phase 
as that in FigEfc) of model 1. The surface in FigOfd) is 
separated from that in FigOJb) by a first-order transition 
as we will see below. 

From the bending energy S2 we can extract informa- 
tion about how smooth the surface is, although S2 is not 
included in the Hamiltonian. The bending energy S2 is 
defined by 

S 2 = ]T (1 - n, • n,) , (7) 

ij 

where is the unit normal vector of the triangle i. The 
symbol ij in J^ij denotes the nearest-neighbor triangles 
i and j, which have a common bond. The pair • n., is 
not defined on the boundary bonds, because we have no 
triangle outside the boundary. 

Figures E{a) and E{b) are plots of the bending energy 
S^/Nb against a obtained in model 1 and model 2, respec- 
tively. Nb is the total number of internal bonds. We find 



(a) Model 2 (a = 10) (b) Model 2 (a = 32) 




(c) Model 2 (a = 36) (d) Model 2 (a = 44) 



Fig. 3. Snapshots of model 2 surfaces of size N = 5167 at 
(a) a = 10 (crumpled), (b) a = 32 (non-smooth), (c) a = 36 
(transient surface between non-smooth and smooth), and (d) 
a = 44 (smooth). The surface in (c) is almost smooth. Figures 
are shown in the same scale. 




20 „ 40 20 40 



Fig. 4. S2/NB vs. a obtained in (a) model 1 and (b) model 2. 
Nb is the total number of internal bonds. The size of surfaces 
are N = 817, N = 1657, N = 3397, are N = 5167. Dashed lines 
denote the phase boundaries of the models on the N = 5167 
surface. 

from Figs^Ja) andEJb) that S^/Nb is obviously discon- 
tinuous at finite a. Thus, we have shown that both mod- 
els undergo a phase transition between the smooth phase 
and the non-smooth phase. As we have seen in the snap- 
shots, both model have three different phases; the smooth, 
the non-smooth, and the crumpled, which is expected to 
emerge at sufficiently small a. 
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One dashed line in the figures denotes the phase bound- 
ary between the smooth phase and the non-smooth phase 
on the N — 5167 surface, and the other one denotes the 
phase boundary between the non-smooth phase and the 
crumpled phase. Since the transition point of the N — 3367 
surface is almost identical to that of the N = 5167 surface, 
the transition seems not to be influenced by the finite-size 
effects in those system size. We note also that the non- 
smooth surfaces do not appear between the smooth phase 
and the crumpled phase on the N = 817 and N = 1657 
surfaces. The reason of this seems due to the finite-size 
effect. 

The order of the transitions between the smooth phase 
and the non-smooth one is considered as first order, be- 
cause S2/NB is discontinuous as shown above. We re- 
call that a phase transition is called a first-order (or dis- 
continuous) one if some physical quantity discontinuously 
changes even when the quantity is not included in the 
Hamiltonian. On the contrary, we can see no transition 
between the crumpled and the non-smooth phases. In fact, 
S^/Nb and some other physical quantities smoothly vary 
at the transition region of cv, just as in the spherical sur- 
face model [T3] . 




model 2 
N=3367 



3 xlO 3 



Ct=32 

ct=40 (non-smooth) 
(smooth) 



sweeps 



1 1.5 xK 

sweeps 



Fig. 5. Variation of S2/NB against MCS obtained on model 
1 surface of size (a) N = 5167, (b) N = 3367, and on model 2 
surface of (c) TV = 5167, (d) 7V = 3367. 



In order to see how is the convergence of the simula- 
tions close to the transition point, we plot the variation of 
Si/Nb against MCS including the thermalization MCS in 
FigsOJa)-Ei;d): FigEJa) is S 2 /N B obtained at a = 22(non- 
smooth) and a = 24(smooth) on the N = 5167 surface of 
model 1, FigOJb) is those obtained at a = 25 (non-smooth) 
and a — 22 (smooth) obtained on the TV = 3367 surface of 
model 1, FigOJc) is those obtained at a = 32 (non-smooth), 



a = 36 (transient), and a = 40 (smooth) on the N = 5167 
surface of model 2, and Fig. Eld) is those obtained at 
a = 32 (non-smooth) and a — 40(smooth) obtained on the 
N = 3367 surface of model 2. We immediately see that 
the total number of MCS after the thermalization in the 
non-smooth phase is relatively smaller than that in the 
smooth phase in the 7V = 5167 surface shown in FigsEJa) 
and|I2c). The reason of this is because the value of Si/Nb 
is stable in the non-smooth phase after the thermaliza- 
tion. On the contrary, relatively large number of MCS 
was performed in the smooth phase close to the transi- 
tion point, because a smooth surface can collapse into a 
non-smooth surface after a long MC run as stated in the 
final paragraph of the previous section. The variation at 
a = 36 (transient) in FigOJc) indicates that the surface 
seems belong to the smooth phase after long MC run, al- 
though the value of Sz/Nb still remains relatively larger 
than that of a = 40 (smooth). 

We found no variation of S2/N B against MCS indicat- 
ing jumps from the smooth phase to the non-smooth phase 
and vice versa. However, the value of S^/^Vb in the non- 
smooth phase is quite distinct from that in the smooth 
phase after the thermalization. Therefore, the gaps (or 
discontinuities) in S2/NB shown in Figs^Ja) andQ{b) are 
considered as a sign of a discontinuous transition. 
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Fig. 6. X 2 vs. a obtained in (a) model 1 and (b) model 2. The 
size of surfaces is jV=817, iV = 1657, iV = 3367, and TV = 5167. 
Dashed lines denote the phase boundaries of the models on the 
iV = 5167 surface. 



In order to see the size of surfaces, we calculate the 
mean square size X 2 defined by 



X v 



(8) 



where X is the center of the surface. 

We plot X 2 against a in Figs||Ja) andHJb) obtained in 
model 1 and model 2, respectively. The dashed lines drawn 
vertically in the figures denote the phase boundaries of 
the iV = 5167 surface. These lines were drawn at the same 
position as the corresponding dashed lines in Figs^Ja) and 
0Jb). We see in FigUJa) that X 2 discontinuously changes 
at finite a in model 1. The transition point a, where X 2 
discontinuously changes in FigOJa), is identical to a where 
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S2/NB discontinuously changes in Fig0{a). The transition 
point a shown in FigOJa) distinguishes the smooth phase 
from the non-smooth phase in model 1. This indicates that 
the surface in the non-smooth phase is almost crumpled. 

On the contrary, these two phases can hardly be distin- 
guished by X 2 in model 2, because X 2 in the non-smooth 
phase is almost comparable to X 2 in the smooth phase 
close to the transition point. This phenomenon is appar- 
ent on the iV> 1657 surfaces and consistent with the fact 
that the surface of the non-smooth phase in model 2 is 
relatively oblong and extended as shown in Fig|^b). 



<X 2 > <X 2 > 




Fig. 7. (a) Log-log plots of (X 2 ) vs. TV (o) obtained at the 
smooth phase close to the transition point of model 1, (b) log- 
log plots of (X 2 ) vs. N obtained in the smooth phase (o) and 
in the non-smooth phase (A) close to the transition point of 
model 2. The straight lines and the dashed line are drawn by 
fitting the data to Eq. ©. 

Figure [7Ja) shows log-log plots of (X 2 ) against N, 
where (X 2 ) are obtained by averaging X 2 over several 
points a close to the transition point with the weight of 
inverse error; a — 24, a = 28 for N = 5167, a = 25, a = 30, 
a = 35 for N = 3367, a = 27, a = 28, a = 30 for N = 1567, 
and a = 37, a = 40 for N = 817. The reason why we need 
such averaging of X 2 is that X 2 changes with a in the 
smooth phase in contrast to the case of spherical surfaces 
30,31 . The straight line is obtained by fitting the largest 
three X 2 to 

X 2 ~ N 2 ' H , (9) 
where H is the Hausdorff dimension. Thus we have 

H = 2.12 ±0.25, (smooth) (model 1). (10) 

The result if = 2. 12(25) in Eq. JUJl is consistent with that 
the surface of model 1 is almost spherical in the smooth 
phase, as we have seen in FigEfc). 

The Hausdorff dimension in the branched polymer phase 
is expected to be H = 2, however, we can not extract in- 
formation about H in the non-smooth phase from our nu- 
merical results of model 1. One of the reason is because 
the non-smooth phase appears only on the surface of size 
N = 3367 and N — 5167, and another reason is because 
X 2 is not always constant against a in the non-smooth 
phase on those surfaces, as we can see in FigUJa). If we 
could have performed very high-statistics simulations on 



larger surfaces, H can be obtained even in the non-smooth 
phase. 

The solid line in Fig[7fb) is drawn by fitting (X 2 ) in 
the smooth phase of model 2; they were obtained by av- 
eraging X 2 close to the transition point; a = 40, a = 44 
for iV = 5167, a = 40, a = 48 for iV = 3367, a = 40, a = 44, 
a = 48 for N = 1567, and a = 55, a = 60 for N = 817. Then, 
we have 

H= 1.73 ±0.17, (smooth). (11) 

The result H= 1.73(17) is considered to represent a prop- 
erty characteristic to the smooth phase as in the case of 
model 1. In fact, H = 1.73(17) in JTTJ is close to H=2, and 
therefore it is consistent with that the surface of model 2 
is almost spherical in the smooth phase, as we have seen 
in FigCUd). 

The dashed line in Fig|7Jb) is drawn by fitting (X 2 ) in 
the non-smooth phase of model 2; they were obtained by 
averaging X 2 close to the transition point; a = 24, a = 28, 
a = 32 for N = 5167, a = 32 for N = 3367, a = 32, a = 36 
for N= 1567. Then we have ff = 1.27±0.09, which implies 
that the non-smooth phase seems not to be the branched 
polymer phase at least in model 2. The non-smooth phase 
seems rather close to an oblong linear surface, because 
H= 1.27 ±0.09 is close to if = 1. 




20 „ 40 20 40 



Fig. 8. (a) Si/N vs. a for model 1, and (b) Si/N vs. a for 
model 2. Si/N has a clear jump in (b). 

Figures E[a) andOJb) show the Gaussian bond poten- 
tial Si/N of model 1 and model 2, respectively, although 
Si is not included in the Hamiltonian of model 2. We 
find from 

Fig0» that the relation Si/N = 3/2, which 
is expected from the scale invariant property of the par- 
tition function of model 1, is not influenced by the first- 
order transition. Whereas in model 2 S±/N discontinu- 
ously changes at the transition point, where S2 has a gap 
as shown in FigQfb). This discontinuity of Si also indi- 
cates that model 2 undergoes a discontinuous transition. 
Therefore, the discontinuous transition of model 2 is char- 
acterized by gaps of £2 and Si, although X 2 seems contin- 
uous at the transition point as we have seen in FigEJb). 

The intrinsic curvature energies S3/N of model 1 and 
model 2 are shown in FigsOJa) and E{b) in a log-linear 
scale. The discontinuity of S3/N is relatively unclear com- 
pared with that of S^/Nb in FigsQJa)BJb). Nevertheless, 
S3 [N also changes discontinuously at the transition point 
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Fig. 9. (a) S 3 /N vs. a for model 1, and (b) S 3 /N vs. a for 
model 2. 



in both models. Moreover, we find also that S3/N in the 
crumpled phase is quite larger than those in the other 
phases. These indicate that that both non-smooth surface 
and crumpled surface are non-flat as an intrisnsic curva- 
ture surface. Thus, the non-smooth surface and the crum- 
pled surface are non-flat both extrinsically and intrinsi- 
cally. 
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Fig. 10. Auto-correlation coefficient A(X 2 ) of (a) model 1 at 
a = 20 (non-smooth), a — 22 (non-smooth), a = 28 (smooth 
phase), and (b) model 2 at a = 32 (non-smooth), a = 36 (tran- 
sient), a — 44 (smooth). 



In order to see the convergence speed of the MC sim- 
ulations in both models, we plot in Figs llOf a) and II Of b) 
the autocorrelation coefficient A(X 2 ) of X 2 defined by 



A(X 



Ei* 2 (7-<)*W) 



H+l 



E<* 2 (r<)r 

n + n x 500, n = 1,2, 



(12) 



The horizontal axes in the figures represent 500 xn (n = 
1,2," -)-MCSs, which is a sampling-sweeps between the 
samples X 2 (n) and X 2 (ri+i). The coefficients A{X 2 ) in 
Figs llOf a) and HOf b) are obtained on the surface of size 
iV = 5617 in model 1 and model 2, respectively. 

We clearly see from Fig llOf a) that the convergence 
speed in the smooth phase (a = 28) is almost identical 
to that in the non-smooth phase (a = 20), and it is also 



identical to that at a = 22, which seems very close to the 
transition point but remains in the non-smooth phase. The 
reason why the convergence speed in the smooth phase is 
comparable to the one in the non-smooth phase is that 
the phase space volume (C R 3 ), where the vertices Xi 
take their values, in the smooth phase is almost identical 
to that in the non-smooth phase. Figure ITUT b) shows the 
coefficient A(X 2 ) of model 2. They were obtained in the 
non-smooth phase (a = 32), in the smooth phase (a = 44), 
and at a transient point (a = 36) between the two phases. 
We find that the convergence speed in the non-smooth 
phase (a = 32) is relatively faster than that in the smooth 
phase in model 2. 



Decorrelation MCS is about 0.3 x 10 7 



1 x 10 7 



the N = 5167 surface close to the transition point in 
both models, because A(X 2 ) - at 0.3 x 10 7 - 1 x 10 7 
in FigsHISfaHToTbl. Therefore, the total number of MCS 
(0.4 x 10 8 ~ lx 10 8 ) after the thermalization, which was 
shown in Figs[5fa) |5fd), is considered to be sufficient, be- 
cause it is very larger than the decorrelation MCS. How- 
ever, we must note that the statistics is still not so high to 
obtain the Hausdorff dimension in the non-smooth phase 
and in the smooth phase as stated before. 

We should mention that phenomena of the critical 
slowing down in A(X 2 ) at the transition point is very 
hard to see in our MC technique. The reason of this is be- 
cause of very low speed of convergence. The change from 
the non-smooth phase to the smooth one appears to be 
almost irreversible. If the vertices are once trapped in a 
small region of R 3 at the transition point, they hardly 
expand to be a smooth surface at least with the simula- 
tion technique in this paper. We obtain no series of MC 
data which includes configurations in the smooth phase 
and those in the non-smooth phase alternately even at 
the transition point. We consider that this irreversibility 
is only due to the simulation technique. These trapping 
phenomena were also seen in the same model on a sphere 
[TU| , which was studied by using the same simulation tech- 
nique. But, the transition between the two different phases 
must eventually be seen after long MC simulations at a 
sufficiently apart from the transition point. 

Finally, we comment on the coordination number of 
surfaces. The maximum coordination number g max is g max = 
22 on the iV = 5167 surface of model 1, and (? max = 19 on 
the N = 5167 surface of model 2. These values of g max 
were almost independent of a at 2 < a < 28 for model 1 
and at 10 < a < 44 for model 2. 



5 Summary and conclusion 

The purpose of this study was to understand how the flu- 
idity of dynamical triangulation influences the first-order 
crumpling transition observed in the model on the fixed 
connectivity disk, which is an open surface. We were also 
interested in that the transition is independent of whether 
the surface is open or closed. 

Two types of fluid surface model have been investi- 
gated on a dynamically triangulated disk with intrinsic 
curvature; one with the Gaussian bond potential and the 
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other with a hard-wall potential in place of the Gaus- 
sian bond potential. Both potentials are defined on all the 
bonds including the boundary bonds. The triangulated 
disk is obtained by dividing the edges of the hexagon into 
L-pieces, for a starting configuration of the MC simula- 
tions. The triangulated disk is characterized by (N, Ne, N't) 
(3L 2 + 3L + 1,9L 2 + 3L,6L 2 ), where N E ,N T are the total 
number of edges and the total number of triangles. In 
order to see whether the finite size effect influences the 
transition or not, we used surfaces of size up to TV = 5762, 
which is relatively larger than those previous studies on 
dynamically triangulated surfaces. The boundary bonds 
are prohibited from flipping, and the remaining internal 
bonds are allowed to flip so that at least one of the two- 
vertices of the bond connects to an internal vertex after 
the flip. 

We have shown that both models undergo a first-order 
phase transition between the smooth phase and the non- 
smooth phase. The transition in model 1 is characterized 
by a discontinuity of 5*2 and that of X 2 , whereas the tran- 
sition in model 2 by a discontinuity of S% and that of S\. 
Those discontinuities are considered to reflect a discontin- 
uous transition. A phase transition is considered to be of 
first-order whenever some physical quantity has a gap or 
a jump. 

The transition occurs independent of whether the Gaus- 
sian bond potential is included in the Hamiltonian or not. 
Surfaces appear to be spherical in the smooth phase in 
both models. Our results indicate that the first-order tran- 
sition is not influenced by the diffusion of vertices over 
dynamically triangulated surfaces if the intrinsic curva- 
ture is included in the Hamiltonian as the curvature term. 
Moreover, the results in this paper together with that in 
indicate that the transition is independent of whether 
the surface is open or closed. 

It is interesting to study the model on the disk with 
grand canonical simulations so that the number of bound- 
ary vertices or bonds can change. 

This work is supported in part by a Grant-in- Aid for 
Scientific Research, No. 15560160. The author (H.K.) ac- 
knowledges the staff of the Technical Support Center at 
Ibaraki National College of Technology for their help on 
computer analysis. 
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